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We find that the frequency spectra of layered phononic and photonic composites admit a universal struc¬ 
ture, independent of the geometry of the periodic-cell and the specific physical properties. We show that this 
representation extends to highly deformable and multi-physical materials of tunable spectra. The structure is 
employed to establish universal properties of the corresponding band-gaps, and to rigorously determine their sta¬ 
tistical and optimal characteristics. Specifically, we investigate the density of the gaps, their maximal width and 
expected value. As a result, rules for tailoring the laminate according to desired spectra properties follow. Our 
representation further facilitates characterizing the tunability of the band-structure of soft and multi-physical 
materials. 


Wave propagation in heterogeneous media has fascinated 
the scientific community for decades. The inhomogeneity 
causes multiple scattering, and, in turn, wave interferences 
that give rise to intriguing phenomena in various fields. Of 
particular interest are the transition of conducting to isolat¬ 
ing behavior of electronic crystals [1], localization of elec¬ 
tromagnetic waves in dielectrics [2], and attenuation of me¬ 
chanical motions [3-8] in elastic media. The significance of 
the latter stems from its central role in numerous applications; 
transducers [9], waveguides [10], vibration filters [11], acous¬ 
tic imaging for medical ultrasound and nondestructive testing 
[12], noise reduction [13] and cloaking [14, 15] are just a few 
examples. The mathematical and physical richness of elastic 
waves in heterogeneous materials emanates from their vecto¬ 
rial nature, and their spatial dependency on additional con¬ 
stituents parameters. 

Layered media have been extensively studied [16-22], in 
virtue of their relative simplicity of fabrication and theoreti¬ 
cal modeling. This letter provides new insights on the relation 
between their geometry, physical properties, and frequency 
spectrum. We find a universal representation for the spec¬ 
trum, independent of the unit-cell geometry and specific con¬ 
stituents properties. This structure reveals a universal property 
of the gaps-density, namely, its invariance under the change of 
various geometric and physical properties. We utilize this rep¬ 
resentation to determine exactly the density of the gaps, and 
their expected and maximal widths. These results are iden¬ 
tified with classes of compositions, hence provide rules for 
tailoring the laminate according to desired spectra character¬ 
istics. Thus far, such calculations would necessitate the trun¬ 
cation of infinite spectra, thereby leading to estimates rather 
than accurate results. 

The canvas upon which the analysis is presented is of 
phononic crystals. The conclusions we draw, however, ex¬ 
tend to additional systems. In virtue of the similarity between 
electromagnetic and elastodynamics wave equations for the 
considered geometry, our insights apply to one-dimensional 
(ID) photonic crystals as well. Stratified piezoelectrics admit 
a similar spectrum too [23]. Our analysis further applies to 
soft non-linear media and multi-physical composites, whose 
physical properties are changed upon application of external 
stimuli. By these means, the frequency spectrum of such ma¬ 
terials is rendered tunable. Our approach facilitates character¬ 



Figure 1. (a) A ID crystal of alternated layers, (b) Attenuating 
waves, associated with frequencies for which \r]\ > 1. 

izing this tunability, as demonstrated in the forthcoming. 

We begin by considering a ID crystal made out of two al¬ 
ternating phases (Fig. la). We denote the phases with 1 and 2, 
and their associated quantities with superscript {p),p = 1 and 
2, respectively. The dispersion relation governing the propa¬ 
gation of waves in the crystal is [16, 20, 24] 

arccos77 = kph^ (1) 

where 

^ = cos^cos^-7Sin^sin^. (2) 

In the above, uj is the frequency, is the Bloch-parameter, 
is the thickness, and h = -b h^‘^\ The parameter 
corresponds to different velocities, depending on the type of 
media and waves considered. Specifically, when the crystal 
is photonic, is the velocity of light in the phase; when 
the crystal is phononic, is the velocity of either trans¬ 
verse waves or longitudinal waves, propagating in a bulk. The 
parameter 7 = ^ ^ quantifies the 

contrast between the constituents impedance, where a = lm 
the photonic case, and a = in the phononic case, 

p^^^ being the mass density. The frequency spectrum is ob¬ 
tained by solving Eq. (1) for values of in the irreducible 

Brihouin zone [25], ^ < kph < tt. Band-gaps correspond to 
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Figure 2. (color online), (a) Representative contours of the disper¬ 
sion relation of real /cs, at 7 = 5, on the unfolded torus, (b) Rep¬ 
resentative domains of D, at exemplary values of 7 = 2,5 and 10, 
denoted by the gray, blue, and red regions, respectively, correspond¬ 
ing to ranges of band-gaps. The red dots designate the intersections 
of D with the ray emanating from the origin, identified with the 

gap. 


ranges of frequencies of attenuating waves, for which |? 7 | > 1 
(Fig. lb). 

Our analysis of the spectrum benefits from a technique 
developed for the study of Schrodinger operators on metric 
graphs [26-28]. The systems we consider necessitate a suit¬ 
able variation of that approach, as follows. Upon defining the 
variables 


(^(p) 


c(p) ’ 


(3) 


one can write 77 as a doubly 27r-periodic function of and 
namely 

V } 7 ) •= cosC^^^ cos^"^^^ — 7sin^^^^ sinC^^^- 

(4) 

Fixing 7, we consider 77 as a function defined on a 2D torus of 
edge length 27r, characterized by the coordinates }• 

The absolute value is invariant under the transformations 
c(^) ^ + TT. This symmetry allows to make a further 

reduction, and fold the torus into a 7r-periodic torus; we denote 
the new torus by T, and denote by B its subdomain where 
|? 7 | > 1 , i.e., where Eq. (1) is solved with imaginary ks- For 
illustration, representative contours of the dispersion relation 
of real ks are plotted in Fig. 2a, on the unfolded torus, at 7 = 
5. Also, representative domains of B, at exemplary values of 
7 = 2,5 and 10, are depicted in Fig. 2b, by the gray, blue, and 
red regions, respectively. 

This geometric representation engenders new insights on 
universal characteristics of the frequency spectrum, as shown 
in the sequel. Towards this end, we interpret Eq. (3) as one 
which defines the following flow on T 


c (w) = w • 


772 22 ) 


mod TT, 


(5) 


where uj has the role of a time-like parameter. Band-gaps are 
identified with values of uo for which the flow (cj) coincides 
with B. We denote the lower and upper curves bounding B by 
Cl and Cu, respectively. We find that 77 = —1 on these curves 


[29] ; plugging this value into Eq. (4) provides the following 
expressions for the curves 


(^(^) = tt — arccos 


cos ±7\/7^ - 1 sin^ 72 

1 + (7^ — 1) sin^ ("B) 


( 6 ) 


where the upper (resp. lower) curve Cu (resp. Ci), cor¬ 
responds to the plus (resp. minus) sign in the numera¬ 
tor. The function 77 on T is invariant under tt -rotation of T 
around its middle, given by the transformation ^ 

(tt — c(i) , TT —and under reflection across the line 
(^( 2 ) = TT — (^( 1 ), by the transformation ^ 

(tt — Each of those transformations leaves 

the domain B invariant and exchanges between its boundary 
curves Ci and Cu- 

The direction of the flow on the torus is given by the ratio 


“ “ hW c(2) • 


(7) 


The irrationality of the ratio, on account of the physical nature 
of the parameters involved, implies that the flow covers the 
torus ergodically, with a uniform measure. Hence, the density 
of the gaps is simply the relative area of B in T, which can be 
calculated via the integral of the closed-form expression 


( 8 ) 

0 


We arrive at a counter-intuitive and peculiar result; at a pre¬ 
scribed 7, the gap-density is independent of the volume frac¬ 
tions of the phases. To shed light on the signiflcance of this 
statement, consider the following example. Compose a lami¬ 
nate of equal volume fractions of materials (1) and (2). Com¬ 
pose a second laminate by introducing an infinitesimal amount 
of (1) into a bulk of (2). The probability that an arbitrary 
frequency pertains to a band-gap is identical in the two lami¬ 
nates. 

The widths of the gaps are investigated next. We recog¬ 
nize that these widths, which we denote by Acc, are related to 
lengths of intervals directed along the flow, which we denote 
by whose endpoints he on Cu and Cu via 

where We associate each length A( with 

the parameters a and b which characterize the line equation of 
its corresponding flow interval 


+ b. ( 10 ) 

These observations, together with the derived expressions for 
the curves Cu and Cu enable determining the width of the gaps, 
and relating it to the physical and geometrical properties of the 
crystal. As the whole spectrum is encapsulated in the torus, 
we are able, in turn, to formulate optimization problems rig¬ 
orously, derived in the sequel. 
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We start with the gap, whose width maximization is of 
practical importance, being the one which is most often re¬ 
alized experimentally [17, 20, 21]. We would like to know: 
given two materials, what is the microstructure which maxi¬ 
mizes the gap? The gap is identified with the flow line 
emanating from the origin. Therefore, we seek the slope a 
which maximizes the right-hand side of Eq. (9), at 6 = 0. The 
problem is interpreted as a search for an optimal at 

fixed . 

To that end, the calculation over the torus serves as an al¬ 
ternative to a calculation of the gap-width for all possible 
compositions, via a partial evaluation of the spectra. Calcu¬ 
lations over the torus become a necessity, however, when the 
whole spectrum needs to be analyzed. Such cases are con¬ 
sidered next, starting with the calculation of the greatest gap 
width. We address two different practical scenarios. The first 
scenario considers a given crystal, with prescribed geomet¬ 
rical and physical properties. The objective is delivered by 
finding the maximal within the segments along flows of 
the prescribed slope a, i.e., maximizing over translations b, 
and plugging it into Eq. (9). We emphasize that this optimiza¬ 
tion is defined solely over D, in an exact manner. Contrarily, a 
direct approach will require the evaluation of an infinite spec¬ 
trum; since in practice it must be truncated, such approach 
provides only an approximation. The second scenario we con¬ 
sider is of a laminate with prescribed constituents, prior to de¬ 
termining their volume fractions. In this case, the maximal 
gap is obtained by maximizing the right-hand side of Eq. (9), 
over all slopes a and translations b. 

We supplement the discussion regarding optimality noting 
that our approach allows for quick calculation of bounds on 
the gaps-width, based on bounding A^. We find that for 7 < 
7 cr, where 7 cr ~ 5.45, the maximal denoted A(^rnax, is 
obtained at the intersection of D with = 7 r/ 2 . Utilizing 
Eq. ( 6 ), it is expressed in closed-form as 


and the variance is calculated by 

b=7T 

b=—a'K 

with E = Ib^-an emphasize again that 

deriving these results using a direct approach involves calcu¬ 
lations over truncations of infinite spectra, thereby leading to 
estimates rather than rigorous results. 

We are ready to introduce novel composites for which our 
analysis applies as well, and begin with laminates comprising 
non-linear soft layers [22, 30]. We consider waves propagat¬ 
ing in such laminates when finitely deformed in a piecewise- 
homogeneous manner. By applying the theory of incremental 
elastic deformations [31], Eq. (1) is recovered as the disper¬ 
sion relation of such “small on large” waves, and therefore the 
resultant spectrum is endowed with the torus representation. 
The physical and geometrical quantities entering Eq. (1) are 
now functions of the finite deformation, which thus renders 
the spectrum tunable. In addition to the validity of the previ¬ 
ous results, our representation further provides a convenient 
tool for characterizing this tunability. This is demonstrated 
next, by way of example. The example considers transverse 
waves superposed on a plane deformation of an incompress¬ 
ible laminate. We denote the principle stretches in each phase 
by and , on account of incompressibility. The lam¬ 
ination direction considered is along the principal axis associ¬ 
ated with such that is related to the thickness before 
the deformation, via The velocity of 

transverse waves propagating in this direction is 


cb) = 



(15) 


ACmax = TT — 2 arccos 



( 11 ) 


When 7 > 7 cr, associated with large impedance contrast, 
ACmax is attained along the diagonal, for which 
Upon substitution of A^^max into Eq. (9), the following bound 
is derived 


Acc < 






=AC„ 


( 12 ) 


Thus, an explicit expression in terms of the crystal properties 
is obtained. So far, to the best of our knowledge, such bounds 
were not accessible. 

The analysis is concluded with expressions for statistical 
characteristics of the spectra, evaluated exactly using our rep¬ 
resentation. Specifically, the expected value of the gap-width 
of a prescribed system is determined via 

( 13 ) 


where the instantaneous stiffness, depends on the under¬ 
lying deformation and the constitutive law of the phase. If a 
neo-Hookean law governs the layers behavior, we have that 

where is the shear modulus in the limit of small strains. 
A simple calculation shows that and the flow direction 
are unaffected by the deformation, namely, 

^ \/Ab)Vb)/pb) ^^{P)/p{P) 

To examine how the domain of D is changed by the defor¬ 
mation, assume without loss of generality that > 

p( 2 )^( 2 ). If ^( 1 ) ^ ^( 2 ) ^( 1 ) then 7 is greater 

(resp. less) than its value before the deformation, and the area 
of B increases (resp. decreases). Together with the observa¬ 
tion made on a and ti^P\ it implies that frequencies are added 
to (resp. removed from) the boundaries of the gaps, without 
changing the mean frequency of each gap. 
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Consider next layers that are characterized by a Gentian 
model [32]. In this case, the instantaneous stiffness is 




Ab)" + A(J’) " - 2"^ 

TW ) ’ 
/ 


(18) 


where the material parameter Jm ^ ( 0 , 00 ] models the rapid 
rise of stress in elastomers when approaching a limiting strain. 
Assume first that = J and = A. In 

such case, the fiow direction and impedance contrast are not 
modified by the deformation; however, the fiow rate, d( /dou, 
is multiplied by 1 — ^A^ + A ^ —2^ jj. Consequently, the 

pertinent frequencies are multiplied by the inverse of this fac¬ 
tor. Observing this factor is a monotonically decreasing func¬ 
tion of A with a range (0,1] [33], implies that the frequency 
spectrum is shifted towards higher frequencies, and its gaps 
are rendered wider, by that factor. When ^ or 
A^^^ 7 ^ A^^\ the fiow direction changes with the deformation, 
and so does the impedance contrast. Thus, the spectrum is 
rendered tunable in a more intricate manner. The exploration 
of the specific way it takes effect depending on the relation 
between and A^^^ is beyond the scope of this example; 
however, we note that following the way the aforementioned 
relations enter Eq. (7), the fiow lines will rotate either clock¬ 
wise or counter-clockwise, and the length of their intersec¬ 
tions with the resultant domain of D will accordingly change. 
The analysis of the resultant Acc is rendered simpler, in virtue 
of the torus universality, as for each 7, the relation between 
and a is unchanged. Therefore, one can calculate the re¬ 
sultant slope a in terms of and A^^\ and employ a fixed 
calculation of AC/Vl + to evaluate Eq. (9). 

We continue with dielectric elastomers (DEs), a class of 
soft active materials. DEs undergo finite strains and change 
their physical properties by application of electric stimuli 
[34, 35]. The work in Ref. [36] has shown that under par¬ 
ticular settings, the propagation of waves in finitely deformed 
DE laminates is also described by Eq. (1). Hence, the torus 
representation holds in this case as well, and subsequently, the 
validity of our previous observations is established. In partic¬ 
ular, the torus representation provides a platform for an in¬ 
vestigation of the effect the electric field has on the frequency 
spectrum. 

We consider next magnetorehological elastomers (MREs) 
[37]. These materials consist of ferromagnetic particles em¬ 
bedded in a rubber-like matrix. Application of magnetic stim¬ 
uli induces magnetic forces and moments on the inclusions. 


This changes the microstructure of the material, and, in turn, 
alters its configuration and stiffness. When the particles are 
disturbed in a chain-like manner, the material admits a lami¬ 
nated structure [38, 39]. We argue that our analysis also ap¬ 
plies to MREs, in light of the similarly with DEs. Eirst, we 
note that under a quasi magneto/electrostatic approximation, 
the governing electric and magnetic fields are differentially 
similar [40] . Second, the mechanical response of both MREs 
and DEs is governed by an elastomeric substance. The only 
difference is the relation between the magnetic load and the re¬ 
sultant stretch, on account of a different magnetic constitutive 
behavior. Therefore, under the same settings, the propagation 
of superposed waves is governed by the same dispersion rela¬ 
tion derived for DEs. The applicability of our analysis thereby 
follows. 

We complete this letter with a summary of our main conclu¬ 
sions, and a glance towards future challenges. We found that 
the frequency spectra of various ID crystals admit a univer¬ 
sal structure, independent of the geometry of their unit-cells 
and specific physical properties. This structure enabled us to 
derive universal properties of the spectra, and rigorously de¬ 
termine the gaps-density, their expected and maximal widths, 
and relate these to particular compositions. We showed that 
our conclusions successfully apply to certain multi-physical 
materials, of tunable spectra. We utilized our framework to 
characterize this tunability. Thus far, the results above were 
either unknown, or determined in an approximate manner. 

Our future objectives are to extend this approach to an¬ 
alyze interesting generalizations of the considered systems. 
One such generalization is a crystal consisting of more than 
two constituents. In this case, the dispersion relation will de¬ 
pend on additional products [41] in the form of Eq. (4), as in¬ 
ferred from a transfer-matrix analysis [36, 41]. Our extended 
analysis will consequently require a torus whose dimension¬ 
ality equals the number of constituents. N-periodic media are 
also systems of major interest. In these composites, a closed- 
form expression for the dispersion relation is not available, 
and series-type solutions are sought [3, 7]. It is an impera¬ 
tive challenge to establish counterparts of the results reported 
herein for such materials. 
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